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ABSTRACT 

We present a kinematic analysis of the globular cluster systems and diffuse stellar light of four 
intermediate luminosity (sub-L*) early-type galaxies in the Virgo cluster based on Gemini/GMOS 
data. Our galaxy sample is fainter (—23.8 < Mr < —22.7) than most previous studies, nearly 
doubling the number of galaxies in this magnitude range that now have GC kinematics. The data 
for the diffuse light extends to 4 R e , and the data for the globular clusters reaches 8-12R e . We find 
that the kinematics in these outer regions are all different despite the fact that these four galaxies 
have similar photometric properties, and are uniformly classified as “fast rotators” from their stellar 
kinematics within 1 R e . The globular cluster systems exhibit a wide range of kinematic morphology. 

The rotation axis and amplitude can change between the inner and outer regions, including a case 
of counter-rotation. This difference shows the importance of wide-field kinematic studies, and shows 
that stellar and GC kinematics can change significantly as one moves beyond the inner regions of 
galaxies. Moreover, the kinematics of the globular cluster systems can differ from that of the stars, 
suggesting that the formation of the two populations are also distinct. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: kinematics and dynamics - 
galaxies: star clusters: general 


1. INTRODUCTION 

Early-type galaxies (ETGs) are thought to be the final 
products of mergers and accretion processes in a hier¬ 
archical universe. This history of collapse and merg¬ 
ing is preserved in the motions of their stars. Stel- 
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lar kinematics can therefore be used to constrain galac¬ 
tic formation and evolution. Absorption line profiles 
measured from integrated spectroscopy can measure 
these bulk stellar motions. However, these observa¬ 
tions are usually restricted to wit hin 1-2 effective radii 
(R < 1-2 R P ) of galaxy centers (Gerhard et al.l 120011: 


Pinknev et al’l 120031 iProctor et al.l 120051 : ISpolaor et al l 
2008f) due to the faintness of the stellar halo. With 


the development of integral field units (IFUs), this 
technique has given us a new two-dimensional view 
of the line -of-sight velocity distribution (LOSVD 


Sanchez et al.1 2012 

: Emsellem et al. 

2014; Bundv et al.l 

20151 Guerou et al. 

201511. iCappellari et al.l (|2011|) ob- 


served 260 ETGs in the nearby universe with the Spec- 
trographic Areal Unit for Research on Optical Nebu¬ 
lae (SAURON) IFU in the ATLAS 30 project. They 
were able to measure the projected angular momentum 
per unit mass, and used this parameter to classify 260 
ETGs as either “fast” or “slow” rotators, finding that 
nearl y all but the most m assive ETGs are fast rota¬ 
tors (jEmsellem et all 1201 111 . They argue for a shift in 
the existing paradigm for ETGs, where galaxies gener¬ 
ally classified morphologically into disc-like SO galaxies 
and spheroidal-like E systems should instead be classified 
kinematically. 

However, like earlier studies, the ATLAS 30 project is 
limited to 1 R e , where relaxation processes make infer¬ 
ring formation histories increasingly difficult, and where 
feedback processes complicate things even further. The 
fraction of the total mass in this region may only be 
~15-30% i Wn et all 120141) . In addition, the outer re¬ 
gions may present more obvious signatures of merging or 
accretion events. There are only a few studies that ex- 
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tend out to ~2-4 R e (Mehlert et al.l 120001: iProctor et al.1 


20091: iWeiimans et al 


200S~ Forestcll &: Gebhardtl [20101: 


Arnold et al.112014 . IProctor et al.l (120091) developed a 

new technique for extracting integrated spectra of the 
field stars, or galaxy diffuse light (GDL), contained in 
multi-slit observations, and found a variety of rotation 
profiles beyond 1 R e for five ETGs. They found that the 
stellar kinematics in their galaxies’ inner regions were not 
necessarily predictive of those in the outer regions. One 
possibility is that dissipation processes in gas-rich major 
mergers produce dramatic ki nematic differences be tween 
the inner and outer regions (jHoffman ct al.l [20101 ). Al¬ 
ternatively, the galaxies could have formed in two sepa¬ 
rate phases, with the inner region forming early by dis¬ 
sipation, and the ou ter region formin g later by accretion 
of smaller galaxies (IQser et al.lI2010T) . Further studying 
stellar kinematics at large radius is clearly an important 
direction to pursue. 

Compared to stellar kinematics, which are limited by 
low surface brightness, observations of globular clusters 
(GCs) can extend out to much larger distances (^10i? e ) 
from the galaxy center. Because GCs are compact and 
luminous, they are easily observed in the halos of galax¬ 
ies. GCs are old, coeval collections of stars that can 
survive a Hubble time, and they preserve information on 
the assembly history of their host galaxies. 

In recent years, large samples of GC vel ocities around 
ETG s have been produced (e.g., M87, Strade r et all 
MTU NGC 1399 ISchuberth et all IMlfl: NGC 5128, 
Peng et al,; 21)1) 111 iWoodlev et all 120101! . For practical 
reasons, much of the work in this field has focused on 
massive early-type galaxies and their rich GC systems. 
ETGs at intermediate mass, however, are much more 
common, and perhaps less complicated than their more 
massive counterparts. For these ETGs, with luminosi¬ 
ties at or fainter than the kn ee of the galaxy luminosity 
function ( Mk = —24.2 mag. iCole et al.ll200lL h — 0.7'). 
there are only a f ew studies (NGC 3115, Arnold et all 
Mil NGC 3379. IPuzia et al.ilMoi IPierce et, all " 


Bergond et al.ll200fil NGC 44 94.IFoster et al. 1120111 NGC 
821, 3377, 1400, 4278, 7454, IPota et all 120131) . Due to 
lower surface brightnesses and sparser GC systems, these 
g alaxies generally lack kinem atic data at large radii. 

iRomanowskv et abl (l2003h and iDouglas et al.l (I2007T ) 
used planetary nebulae as kinematic tracers and found 
some intermediate luminosity ETGs (NGC 3379, NGC 
821 and NGC 4494) with low mass-to-light ratios , which 
appeared to conflict with the dark matter fractions pre¬ 
dicted in cosmological simulations. Compared to mas¬ 
sive ETGs, intermediate luminosity ETGs may also have 
more sporadic merger histories with different mass ra¬ 
tios, potentially preserving the wide range of assembly 
processes at work, and making them important and in¬ 
teresting targets of study. 

In this paper, we study the GC and stellar kinematics 
of four intermediate luminosity ETGs in the Virgo cluster 
out to large radii. We compare GC and stellar kinematic 
trends (rotation amplitude and position angle) with each 
other, and at different radii. In our next paper (Li et 
al., 2015b, in prep), we will study the mass distribution 
and dark matter content of these intermediate luminosity 
ETGs. 


2. SAMPLE, OBSERVATIONS, AND DATA 


2.1. Sample 

We selected four intermediate-luminosity ETGs from 
the A CS Virgo Cluster Survey (ACSVCS; ICote et abl 
12004 . a homogeneous Hubble Space Telescope survey 
of 100 ETGs in the nearby Virgo cluster of galaxies us¬ 
ing th e Advanced Camera for Surveys fACS: lFbrd et al.1 
119981) . These deep images, with their high spatial res¬ 
olution^ allow forexcellent selection of globular clusters 
(iJordan et al.ll2009tl . This is crucial for our sample galax¬ 
ies, as they do not have very rich GC systems, and con¬ 
tamination from foreground and background sources can 
significantly reduce the spectroscopic yield. 

The four galaxies in our study are VCC 1231, 2000, 
1062, 685 (NGC 4473, 4660, 4442, and 4350, respec¬ 
tively, see Table [T|). They span the interesting L* to 
sub-A* luminosity range (—20.5 < Mb < —19), making 
them fainter than galaxies targeted in typical GC kine¬ 
matic studies. For its luminosity range, this sample now 
doubles the number of galaxies with GC kinematics. Ac¬ 
cording to NED, two have morphological type E5, and 
the other two are classified as SO. All g alaxies are clas¬ 
sified as “fast rotators” by ATLAS 30 . (lEmsellem et all 

EHID 

2.2. Observation 

We observed these galaxies with the Gemini Multi- 
Object Spectrographs CGMOS. iHook et all 12004 . twin 
instruments on the Gemini North and Gemini South tele¬ 
scopes. Our target galaxies have sizes ( R e ~ 10-18") 
that fit well within the GMOS field-of-view (5.5 arernin 2 ), 
providing coverage out to 10-16A e . Each galaxy con¬ 
tained ~ 50 targetable GCs with V < 23 mag. VCC 
1231, VCC 1062 and VCC 2000 data were taken with 
GMOS-South, whereas data for VCC 685 was taken with 
GMOS-North. We observed three masks for each galaxy 
to overcome the issue of slit-crowding at the galaxy cen¬ 
ters. VCC 685 was only observed with one mask due to 
it being lower priority in the queue. Tabic [5] summarizes 
the observing information. Because there are two gaps 
in the GMOS CCD focal plane, we observed with two 
wavelength centers, 5080A and 5120A, dithering expo¬ 
sures between the two. For each central wavelength, we 
took two exposures. We used the B600_G5323(GMOS 
North) and B600_5303(GMOS South) grisms with 1"0 
wide slits, giving a dispersion of 0.5 and 0.45 A pixel -1 , 
respectivel y. Both grisms gi ve a spectral resolution of 
R = 1688. HPnzia et abllMl 3?) 

2.3. Data Reduction 

The data reduction was performed by using the Gem- 
ini/GMOS IRAF package (Version 1.11). (1) The cos¬ 
mic ray rejection was done using the GSCRSPEC script, 
which was originally written by Bryan Miller. GSCR¬ 
SPEC is an MEF wrapper for the spectroscopic ver¬ 
sion of the Laplacia n Cosmic Ray Identification routine 
(Ivan Dokkuml IMoIf ) : (2) The bias subtraction and flat- 
fielding were done by the GSDEDUCE; (3) We used 
GSWAVELENGTH to establish the wavelength calibra¬ 
tion from the CuAr arc frames. Typical wavelength so¬ 
lution residuals were ~ 0.1A-0.3A; (4) The wavelength 
calibration was applied to the object frames by using 
GSTRANSFORM; (5) We used GSSKYSUB to do the 
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Table 1 

Sample properties 


ID 

(1) 

Type 

(2) 

R.A.(J2000) 

(deg) 

(3) 

Decl.(J2000) 

(deg) 

(4) 

Vsys 

(km s — *) 

(5) 

Re 

(kpc) 

(6) 

D 

(Mpc) 

(7) 

M k 

(mag) 

(8) 

PA k 

(degree) 

(9) 

(b/a) K 

(10) 

VCC1231 (NGC4473) 

E5 

187.453659 

+13.429320 

2260 

1.2 

15.2 

-23.77 

275 

0.54 

VCC2000 (NGC4660) 

E5 

191.133209 

+11.190533 

1087 

0.8 

15.0 

-22.69 

280 

0.5 

VCC1062 (NGC4442) 

SB0 

187.016220 

+09.803620 

547 

1.3 

15.3 

-23.63 

82 

0.51 

VCC0685 (NGC4350) 

SAO 

185.990891 

+16.693356 

1210 

0.9 

15.4 

-23.13 

30 

0.3 


Note. — Hubble types (2) are from the NED database. The position of galaxy center, RA (3) and DEC (4), the galaxy 
systemic velocity (5) and the K band absolute magnitud e (8) are all from A TLAS 3D project jCappellari et ajjj20 LI[), where 
(8) ar e derived from the 2MASS e xtended source catalog dJarrett et al.|[2000D . The effective radii (6) are from Ferrarese et al.| 
(120061). The dist ances (7) are from [Blakeslee et al.l (j2009l) . except VCC 685 which is from the GC luminosity function e stimates 
by I Jordan et aTI d2007D . The photometric position angle (9) and axis ratio (10) are from 2MASS dSkrutskie et ah||2006t) . 


Table 2 

Summary of Gemini/MOS Observation 


Mask 

R.A.(J2000) 

(deg) 

Decl.(J2000) 

(deg) 

Exp. Time 
(sec) 

Slits 

Night 

GS2008AQ008-01 

191.134109 

+11.176211 

3040 

22 

2008-04-06 

GS2008AQ008-02 

191.134109 

+11.176211 

3040 

23 

2008-04-07 

GS2008AQ008-03 

191.134109 

+11.176211 

3040 

20 

2008-04-08 

GS2008AQ008-04 

187.453537 

+13.427367 

3040 

30 

2008-03-14 

GS2008AQ008-05 

187.453537 

+13.427367 

3040 

25 

2008-04-01 

GS2008AQ008-06 

187.453537 

+13.427367 

3040 

27 

2008-04-05 

GS2008AQ008-07 

187.024216 

+09.806433 

2700 

21 

2008-03-02 

GS2008AQ008-08 

187.024216 

+09.806433 

2700 

21 

2008-03-09 

GS2008AQ008-09 

187.024216 

+09.806433 

2700 

21 

2008-03-14 

GN2008AQ073-01 

186.002594 

+16.686339 

2820 

24 

2008-03-31 


sky-subtracted interactively and GSEXTRACT to ex¬ 
tract ID spectra from each object frame; (6) Finally, 
we used the SCOMBINE to median combine extracted 
spectra from the same slits. 

After extracting the target spectrum from each of the 
mask slitlets, we get a “background” spectrum for each 
slitlets. For slits closer to the galaxy center, this “back¬ 
ground” spectrum contains not only sky but also useful 
amounts of host GDL. If the slit is far from the galaxy 
center, the spectrum should be dominated by “pure sky”. 
In order to recover information from the slits with sig¬ 
nificant amounts of galaxy light, we carefully choose slits 
at the greatest distance (~7 R e ) from the galaxy center 
and combine their spectra to produce a pure sky spec¬ 
trum for non-local sky subtraction. We do not use slits 
located near the edges of the mask where the vignetting 
effects are serious. After subtracting the sky spectrum 
from the remaining “background” spectra we are left 
with the GDL spectra. Figure [1] shows GC and GDL 
spectra with different signal-to-noise ratio. Our spectra 
typically cover a wavelength range from 3500A to 6500A, 
but the wavelength range depends on the slit position on 
the mask and some spectra can reach wavelengths as red 
as 7200A. 

2.4. Radial Velocity Determination 

We use the FXCOR task in IRAF to measure the 
GC and GDL velocities. The template we use for 
FXCOR is a K-type M31 globular cluster spectrum 
from a list of Smithsonian Astrophysical Observatory 
(SAO) spectral templates. For the GDL, we also use 
the penali zed pixel-fitting method ( pPXF) as imple¬ 
mented bv lCappellari fe Emselleml (|2004T ) to measure ra- 



4200 4400 4600 4800 5000 5200 5400 5600 
Wavelength(A) 


Figure 1. Reduced GC and GDL spectra. The red lines are the 
best-fitt PPXF templates. The signal-to-noise (S/N) is calculated 
from 4500A to 4750A. 

dial velocities and velocity dispersions. For pPXF, we 
use the MILES li brary of ste llar spectra as templates 
(jSanchez-Blazauez et all 120061) . The radial velocities 
measured using FXCOR and pPXF are very consistent, 
so we choose to use the FXCOR velocities. The veloc¬ 
ity uncertainties are estimated by FXCOR and scaled 
to those determined with repeated measurements of sci¬ 
ence targets. For velocity dispersion, we use the region 
of 4850A to 5350A which contains the H/3, Mg b and 
Fe5270 spectroscopic features. To be conservative, we 
only measure the velocity dispersion within 3 R e . The 
signal-to-noise of the spectra are typically too low to 
measure accurately the velocity dispersion in the outer 
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Figure 2. A comparison of our GDL radial velocity and velocity 
dispersion measurements with those from the SAURON survey. 
The solid line is a 1:1 relation. The mean offset is 16 km/s and 
4 km/s for the velocity and velocity dispersion, respectively. The 
two surveys show good agreement. 

regions. 

We divide our GCs into red and blue popul ations ac¬ 
cordin g to the (g — z)o colour distributions in iPeng et all 
(120061 ). The colors of GCs are from the A CSVCS and the 
Next Generation Virgo Survey (NGVS, iFerrarese et al.l 
120121) . The NGVS data are used for objects that fall 
outside the ACS field of view. The transformation that 
we use between the ACSVCS and NGVS filters are: 

9acs 9ngvs = 0.12 A 0.06 (1) 

z acs zngvs = 0.06 A 0.10 (2) 

The number of spectroscopically confirmed blue and red 
GCs in our galaxy sample is listed in Table [3] Most of 
the GCs belonging to VCC 685 and VCC 2000 are blue. 

In addition to the Gemini/GMOS spectra of the GCs 
and GDL, all four of our sample g alaxies have pub li shed 
IFU spectroscopy from SAURON (lEmsellem et al.l[2004 
ICappellari et alil2011[) . Since we have some GDL slits 
in the SAURON regions, we compare our measured ra¬ 
dial velocity and velocity dispersion with those obtained 
with SAURON (Figure 0- The velocities and velocity 
dispersions of SAURON data are from a 2D interpola¬ 
tion of their radial velocity and velocity dispersion held. 
The solid line is a 1:1 relation. We found generally good 


Figure 3. Velocity comparison between the GDL and the target 
observed in the same slits. We see the data point is random, which 
suggest there is no correlation between the velocity measured for 
the target and the GDL. 

agreement between these two data sets. SAURON ve¬ 
locities are a little bit larger than our GDL velocities. 
The offset is 16 km/s and 4 km/s for the velocity and 
velocity dispersion, respectively, between our GDL data 
and SAURON data. To facilitate comparison to the 
SAURON data, we add a 16 km/s velocity offset to our 
measured velocities. We do not change our velocity dis¬ 
persion measurements, as 4 km/s is well within the esti¬ 
mated uncertainties. 

One possible problem of obtaining stellar kinematics 
with the GDL technique is that the “background” spec¬ 
trum may be contaminated by residual GC light. To 
check whether this is a problem, we compare the veloc¬ 
ity of the GC and GDL from the same VCC 1231 slits. 
In Figure |3l we show that there is no correlation between 
these two velocities, suggesting that the effect of the GC 
light on the GDL is negligible. 

We present the full table of GC and GDL positions, ve¬ 
locities, velocity dispersions (GDL only), and uncertain¬ 
ties in Tables [8] and [9l The radial velocity distribution is 
shown in Figure 01 

3. KINEMATICS ANALYSIS 
3.1. Two-dimensional Velocity Field 

For each galaxy, we show the two-dimensional velocity 
fields as derived from GCs (e.g., Left panel of Figure 0]). 
We use a Gaussian kernel function to smooth the velocity 
field. At each point, the value of the velocity held is 
estimated by: 


V' ( x,y') 


Ej v jUj 

Ej “j 


1 -D 2 


( 3 ) 

( 4 ) 


where ujj depends on D, which is the distance from j-th 
point to (x , y ), and cr which is the width of the Gaussian 
kernel. The Gauss ian kernel can have a fixed bandwidth 
( Peng et al.ll2004bD o r vary as a function of local density 
( Coccato et al.ll2009l ). The Gaussian kernel we used has 
a fixed width of o — 1 kpc, which is close to 1 R e for 
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Figure 4. The GC and GDL velocity distribution as function of radius. The blue, red circle, and green stars represent blue GCs, red 
GCs, and GDL, respectively. The galaxy systemic velocity is plotted as a dashed line. 


our galaxy sample. We note that our final result is not 
strongly sensitive to the exact kernel width, so we choose 
1 kpc as a reasonable compromise. We also found that 
the fixed kernel did not lead to information loss in the 
outer region of galaxy. In the very inner regions, where 
it might matter, we typically have few data points due 
to slit crowding. 

To reduce the noise due to the finite spatial sampling 
of our GCs and GDL, we assume the galaxy is point- 
symmetric in phase space, a valid assumption for tri- 
axial potentials. Each point (x,y,v) thus has a mirror 
counterpart (—x, —y, —v), where the x and y is the po¬ 
sition of the GC/GDL on the sky, v is the radial ve¬ 
locity of GC/GDL minus the galaxy’s systemic veloc¬ 
ity. T hi s technique has bee n adopted in I Arnaboldi et all 
(|1998f k iPeng et all (l2004b| k iCoccato et al.1 (|2009h . For 
the GCs, we plot the smoothed point-symmetric two- 
dimensional velocity fields for our galaxy sample (e.g., 
left panel of Figure [5]). For the GDL, we only plot the 
point-symmetric two-dimensional velocity fields. In or¬ 
der to compare our GDL data to the SAUR ON data 
(lEmsellem et all l2004i : iCappellari et al.ll2011t h we plot 
the SAURON data in the central regions (e.g., right 
panel of Figure [5]). The raw radial velocity fields with¬ 
out smoothing, and without the assumption of point- 
symmetry (and hence without doubling the data) are 
shown in the Appendix. The point-symmetric velocity 
fields are only used for illustration purposes; all quanti¬ 
tative results in the following sections use the original, 
non-symmetrized data. 


3.2. Rotation Amplitude and Rotation Axis 

One of the main kinematic characteristics we wish to 
explore for halo systems is the degree of rotational sup¬ 
port, if any. To obtain the rotation amplitude and axis 
for the GC systems and stars, we use meth o ds which have 
been used previously ((Foster et all 120111 : iStrader et al.1 
120111 [Pota et alJl201.il) . For this kinematic analysis, we 
use the original data without point-symmetric folding or 
smoothing. For each bin in radius, if we assume the ra¬ 
dial velocity distribution is Gaussian, then the equivalent 
X 2 is: 


X 


2 


E 


(Vi Umod) 

CT p + (^ v i) 


+ Hal + Av?) 


(5) 


Umod — riQ T 


v„ 


/-^ | ( tan(PAj —PA kin ) y 


( 6 ) 


where Vi is the radial velocity, Av is the velocity error, 
tip is the velocity dispersion, V ro t is the rotation ampli¬ 
tude, vo is the galaxy systemic velocity. PA;, PAkin, qkin 
is the position angle of each data point, position angle of 
the kinematic major axis, and the kinematic axis ratio, 
respectively. For the GCs, we minimize the Equation 3. 
For the GDL, the Xo is defined as follows: 


Xo = E 

i 


'Cmod) 

(Ar> 2 ) 


+ ln(Au 2 ) 


(7) 


For both the GCs and GDL, we assume qki n is equal to 
photometric axis ratio of the galaxy light since qkin is dif- 

































6 


ficult to constrain with our data. We also fit with qkin = l 
(round), and the results do not show any significant dif¬ 
ferences. We use Monte Carlo simulations to estimate 
the errors of the fitted parameters. At each data point, 
we generate the velocity with the Gaussian distribution 
defined by the intrinsic dispersion and observed velocity 
error. After 1000 Monte Carlo simulations, we used the 
68 percent confidence intervals as the lcr uncertainties of 
the fitted parame ters. _ 

As discussed in IStrader et al.1 (j2011l h allowing the po¬ 
sition angle of the rotation axis to be a free parameter 
creates a tendency to overestimate the rotation ampli¬ 
tude. We us e Monte Ca r lo sim ulations, similar to those 
described in iPota et al.l (120131b to estimate the magni¬ 
tude of this bias. The bias depends on (V rot /cr), bin size 
and the azimuthal distribution of data points. When 
(V r ot/c) and bin size are small, the bias will increase. 
Sometimes, the estimated bias is comparable to or larger 
than the measured v rot (i.e., for VCC 685). We list 
the estimated bias in Table [3] However, in both cases 
where the bias is comparable to the rotation amplitude, 
the uncertainty in the rotation amplitude is also large. 
Thus, our conclusions are based on the rotation ampli¬ 
tudes without bias correction. 

To compare the SAURON data, GDL, and GCs, we 
carry out rolling fits with the raw data binned in elliptical 
annuli. The number of GCs in each bin depend on the 
sampling of the GC system. There are 10-20 GCs and 
7 10 GDL slits in one bin. We move the bins one data 
point at a time. The galaxy systemic velocity is fixed 
to the value in Table CD The fitting results for the GC 
systems are listed in Table O HI El 

4. RESULTS 
4.1. VCC 1231 

The 2D velocity field of VCC 1231 is shown in Fig¬ 
ure [5] The GC and stellar velocity fields are shown in 
the left and right panels, respectively. The velocity field 
of the GCs does not show any obvious ordered motion 
or significant rotation. The central stellar kinematics as 
shown by the SAURON data, however, show clear rota¬ 
tion, with a peak amplitude of V ro t ~ 60 km s -1 . Our 
GDL data goes out farther in radius than the SAURON 
data and is consistent with the SAURON data in the 
overlap region. In the out er regions, the rota tion am¬ 
plitude declines. Recently, iFoster et al.1 (2013) also use 
a similar method to get GDL kinematics in VCC 1231, 
and their results show both minor and major axis rota¬ 
tion in the outer regions. They argue tha t it is tell-tale 
sign of triaxiality (e.g. JPeng et al.ll2004bf) . Although we 
do not see this kinematic structure in our GDL data, it 
possibly due to the fact that most of the slits where we 
can derive stellar kinematics are along the major axis, so 
it is difficult to see any major axis rotation. We also lack 
the data to confirm the minor axis rotation out to 5 kpc. 

In Figure O we plot the rotation amplitude, V ro t, ve¬ 
locity dispersion, a and rotation axis, PAkin, profiles with 
radius. We can see that a and the kinematic position 
angle of the stars in the outer regions (GDL, green solid 
line or triangle) is generally consistent with that mea¬ 
sured by SAURON in the inner regions (red stars). In 
the region where they overlap, V ro t and PAkin of the 
stars in the outer regions are a little bit offset from the 


SAURON data. One possibility is our GDL kinematics 
are not determined from a uniform spatial distribution, 
which can lead to fitting bias. Moreover, the radial bin 
size for the SAURON profile is much smaller than for our 
GDL profile. The V ro t of GDL data is in agreement with 
SAURON data in 2er confidence. 

The GDL rotation amplitude is smaller than that in 
the inner regions. It decreases with radius and V ro t ~ 
20 km s” 1 at 3 kpc. The PAkin is ~ 120° and stays 
the same out to 3 kpc. Similar results were also found in 
(iFoster et al.ll2013l their Figure 3). We plot the GDL ve¬ 
locity dispersion profile in the middle panel of Figure [6] 
Each green triangle represents the measured velocity dis¬ 
persion from one slit. We derived the SAURON velocity 
dispersion profile (red stars) by taking the mean velocity 
dispersion within an elliptical annulus on the 2D veloc¬ 
ity dispersion map. These two measurements are not 
entirely equivalent, making GDL profile noisier than the 
SAURON velocity dispersion profile. However, on aver¬ 
age, these two datas show a good agreement at the same 
position (see Figure [2]). 

The V ro t of red and blue GCs are similar and both 
nominally larger than that of the GDL, but the error is 
large and the low number of GCs (10 < TV < 20) in each 
bin may lead to the V ro t being overestimated. After the 
bias correction, the rotation velocity of red and blue GC 
is 3±115kms _1 and 52±43kms _1 . The red GC rotation 
is consistent with zero, within the uncertainties. 

Although the rotation amplitudes for the two GC sub¬ 
populations are similar, the best fit rotation axes are 
very different, which is intriguing. The red GCs appear 
to rotate around the photometric minor axis, which is 
consistent with the stars in the inner regions, while the 
blue GCs appear to rotate around an axis somewhere 
between the photometric major and minor axes. These 
very different PAs lead to the whole GC system showing 
no significant rotation (see Figured]). The best-fitting ro¬ 
tation amplitude for the full GC system is 35 ±40 km s -1 
with PAkin = 216 ± 78°. To test the significance of this 
rotation signal, we use the method of “scrambling” the 
properties of the observed GCs and refitting the kinemat¬ 
ics. This method had bee n app l ied i n t esting the rotation 
of the M87 UCD system (IZhang et al.ll2015 l and the GC 
systems of low-mass galaxies. We fix the GC position 
angles and randomly shuffle the velocity and velocity er¬ 
ror of each GC. We repeated this exercise 1000 times and 
each time we apply the same kinematics fitting as we did 
for our original data. We find that 88% and 33% of the 
simulations for the red and blue GCs, respectively, have 
rotation amplitudes at least as high as what we fit for 
our observed data. We conclude that the significance of 
these rotation amplitudes is not high, but future obser¬ 
vations to obtain more GC velocities may shed further 
light on this issue. 

In Figure [7] we also plot the GDL best-fitting rotation 
curve. The GDL best-fitting amplitude is 41 ± 6 km s -1 
with PAkin = 129 ± 7°, which is smaller than the 
SAURON data in the inner region. 

4.2. VCC 2000 

The right panel of Figure [5] shows the 2D velocity field 
of stars in VCC 2000. Both the SAURON and GMOS 
data show obvious rotation with the same PAkin, al¬ 
though the rotation velocity in the inner regions seems 
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Table 3 

All GC kinematics fitting results 


ID 

(1) 

Nblue 

(2) 

N r ed 

(3) 

V m 

(km s -1 ) 

(4) 

Vrot 

(km s -1 ) 
(5) 

Bias 

(km s _1 ) 
(6) 

a 

(km s -1 ) 

(7) 

V rot / O’ 

(8) 

PA ki n 

(degree) 

(9) 

VCC1231 

39 

12 

2268 

35±40 

31 

171±24 

0.2±0.2 

216=1=78 

VCC2000 

42 

1 

1062 

100±42 

5 

128±20 

0.8±0.4 

198=1=32 

VCC1062 

28 

12 

484 

153±28 

8 

108±18 

1.4±0.4 

92=1=12 

VCC685 

15 

0 

1244 

66±120 

106 

195±45 

0.3=1=0.6 

92=1=91 


Note. — GC kinematics results for our galaxy sample. The number of blue GCs, red GCs, the mean 
velocity, rotation velocity (amplitude), bias of the rotation velocity, velocity dispersion, dominance 
parameter (V ro t/cr), and kinematic position angle (PA) are shown in columns (2) to (9), respectively. 
See Section 1 3.2 1 for details on the estimation of these parameters. 


Table 4 

Blue and Red GC kinematics fitting results 


ID 

(1) 

Vrot,B 

(km s -1 ) 
(2) 

V rot ,R 

(km s -1 ) 

(3) 

Biasg 
(km s -1 ) 

(4) 

BiasR 
(km s -1 ) 

(5) 

(km s -1 ) 
(6) 

(km s -1 ) 

(7) 

PAkin.B 

(degree) 

(8) 

PAkin.R 

(degree) 

(9) 

VCC1231 

73=1=43 

83=1=115 

21 

80 

157=1=24 

183=1=51 

235=1=48 

96=1=87 

VCC2000 

97=1=44 

- 

7 

- 

127=1=20 

- 

196=1=30 

- 

VCC1062 

145=1=37 

183=1=57 

9 

20 

108=1=21 

104=1=30 

88=1=16 

100=1=24 


Note. — Blue (B) and red (R) GC kinematics results for our galaxy sample. 


Table 5 

Bright and Faint GC kinematics fitting results 


ID 

(1) 

V rot ,b 

(km s -1 ' 
(2) 

V rot,f 

1 (km s _1 ) 

(3) 

Biasb 
(km s -1 ) 
(4) 

Biasf 
(km s -1 ) 

(5) 

Cb 

(km s _1 ) 
(6) 

CTf 

(km s -1 ) 

(7) 

PAkin,b 

(degree) 

(8) 

PAkin.f 

(degree) 

(9) 

9 

(magnitude) 

(10) 

VCC1231 

118=1=62 

28=1=56 

17 

54 

171±31 

162±30 

224±45 

94±98 

22.274 

VCC2000 

119=1=64 

99 ±46 

20 

16 

138±31 

119 ± 27 

204±45 

182 ± 31 

22.250 

VCC1062 

139=1=29 

170±50 

8 

13 

77±17 

134±28 

88±16 

94±20 

22.458 

Note. — 

Col. (2-9): 

Bright (b) and 

faint (f) GC kinematics results for 

our galaxy 

sample. Col(10): The 

divided g band 


magnitude for the bright and faint group. 


VCC1231 VCC1231 



ARA(kpc) ARA(kpc) 

Figure 5. The 2D velocity fields of our sample galaxies. Left column : Smoothed GC velocity fields. Squares represent the location of the 
GCs, and the triangles represent the point-symmetric counterparts. The color represents the value of the smoothed velocity field at that 
point, with the color scale on the right hand side of each panel. The ellipses represent 1 and 3 effective radii, the values of which are listed 
in Tabled Right column : Stellar velocity field as derived from the GDL data without smoothing. The circles represent the location of 
our slits, and the diamonds represent the point-symmetric counterparts. The SAURON data are plotted in a contiguous rectangle at the 
galaxy center. Note that the spatial and velocity scales differ in each plot in order to best display the dynamic range. In all plots, North 
is up, and East is to the left. 





















R/R. 
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Figure 6. The kinematic fitting results as a function of radius in VCC 1231. The top panel shows the rotation amplitude. The middle 
panel shows the velocity dispersion. The bottom panel shows the kinematic position angle . The red stars represent the fitting results for 
the SAURON data. The lines represent the fitting results from the multislit galaxy diffuse light data (green line), red and blue GCs (red 
and blue), with the shaded regions showing the lcr uncertainty of the fit. The green triangles represent the velocity dispersion measured 
in each individual slit within 3 R e - The photometric major and minor axes are represented with the dotted and dashed line in the bottom 
panel, respectively. The GDL measurements have a rotation axis similar to that of the SAURON data. However, the rotation amplitude 
(top panel) is generally lower than for the SAURON data. Red and blue GCs have the same rotation velocity but rotate around the 
difference axes (although with low significance). The velocity dispersion of the GDL data is consistent with that of the SAURON data. 
The velocity dispersion profile of blue and red GC stays flat out to large distance. 
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VCC 1231 



Figure 7. Radial velocities versus position angle in VCC 1231. 
The black circles represent the GC data and pink triangles rep¬ 
resent the GDL data. The black and pink solid lines represent 
the best-fitting rotation curve for the GCs and GDL, respectively. 
GCs show no rotation or very small rotation. The GDL shows 
small rotation. The curves are not perfect sinusoids in PA because 
we include the ellipticity of the system in the fit. 

larger than in the outer regions. The globular clusters 
show rotation in the inner region < 3 R e and outer region. 
In the middle region, however, the rotation velocity dips 
to ~ 40km s -1 . There are 42 blue GCs and only 1 red 
GC in VCC 2000, so we only plot the blue GCs in Fig¬ 
ure El 

Although there is little overlap between the SAURON 
and GMOS data for the stars, the stellar rotation axis 
(PAkin) is consistent across both data sets, with both 
showing rotation about the photometric minor axis. 
However, the amplitude of the rotation appears to de¬ 
cline with radius. The rotation amplitude of the GC 
system is consistent with that of the stars in their region 
of overlap, but the axis of rotation is off by 90 degrees, 
with the GCs rotating about the photometric major axis. 
The GCs and the stars appear kinematically decoupled 
in this galaxy. The velocity dispersion of the GCs mono- 
tonically declines. 

The best fitting PA versus V ra d curve is shown in Fig- 
ureQjl] The GC best-fitting amplitude is 100±42 km s -1 
with PAkin = 198±32°. The GDL best-fitting amplitude 
is 76 ±6 km s^ 1 with PAkin = 275 ±4°. This again shows 
that the GCs and stars both show rotation, but with the 
different position angles for their rotation axes. The ro¬ 
tation velocity of the GDL is smaller than that derived 
from the SAURON data in the inner regions. 

4.3. VCC 1062 

The 2D velocity fields of VCC 1062 are shown in Fig¬ 
ure mi In this galaxy, both the stars and GCs show 
strong rotation about the minor axis at all radii. The 
blue and red GCs also exhibit similar kinematics. In 
Figure fl2l we can see that the rotation amplitude of the 
stars and GCs at large radii (R > 2 kpc) stays flat at 
around 200 km s -1 . The GCs also have a rotation axis 
similar to that fit to the stellar SAURON data in the 
inner region. 

The best fitting rotation curve of VCC 1062 is shown 
in Figure [13] We can see that the stars (GDL) and GC 
both show a significant rotation with the same PAkin- 


The GCs best-fitting amplitude is 153 ± 28 km s -1 with 
PAkin = 92 ± 12°. The GDL best-fitting amplitude is 
177 ± 3 km s -1 with PAkin = 84 ± 1°. The peak rota¬ 
tion is essentially along the photometric major axis (i.e., 
rotating about the minor axis). The rotation velocity of 
the GDL is larger than that for the SAURON data in the 
inner region. 

4.4. VCC 685 

VCC 685 is also classified as a fast rotator in the 
ATLAS 3D project. The stars clearly rotate ab out the 
photometric minor axis (iKrainovic et all 12011 1. Our 
GMOS GDL measurements are consistent with the 
SAURON data, and we can clearly see rotation in the 
2D velocity field of the stars out to 3 R e (Figure [DO right 
panel, Figure H5l top panel). In these inner regions, the 
GCs also show strong rotation about the photometric 
minor axis, but there is a remarkable shift at R = 3-4 
kpc. At these radii, the rotation axis of the GCs flips and 
exhibits counter-rotation. To better quantify this, we di¬ 
vide the GCs into two groups: 9 GCs in the inner region 
(< 3.5 kpc), and 6 GCs in the outer region ( > 3.5 kpc). 
The best fitting rotation curve for the inner and outer GC 
samples are shown in Figure [16] The purple and green 
circles represent the inner and outer region GCs, respec¬ 
tively. The inner region GCs’ best-fitting amplitude is 
282 it 129 km s -1 with PAkin = 64 ± 38°. The outer re¬ 
gion GCs best-fitting amplitude is 245 ± 85 km s -1 with 
PAki„ = 224 ± 13°. 

Given that we only have 15 GCs in this galaxy, we 
wanted to see if this situation could occur randomly. To 
test the significance of the counter rotation, we used 
the scrambling method discussed earlier where we fix 
the position angles of the GCs and randomly scram¬ 
ble the velocities and velocity errors. After 100,000 tri¬ 
als, we find that only 520 have kinematics similar to 
what we observe (V ro t,inner A 282, V ro t,outer : 245, and 
121° < A PAkin < 239°). With only 0.52% of trials ran¬ 
domly producing the observed kinematics, we conclude 
that this detection of counter-rotation in the GCs is sig¬ 
nificant. 

In Figure Uni we show that the best-fitting amplitude 
of GDL is 174 ± 7 km s -1 with the PAkm is 28 ± 5°. This 
is consistent with the results from fitting the SAURON 
data. These stellar kinematics are roughly (but not ex¬ 
actly) consistent with the kinematics of the inner GCs. 

5. DISCUSSION 

5.1. The Kinematics of GC Systems and Their Hosts 

The four galaxies in our sample are similar in many 
respects. All are classified as “fast rotators” by the 
ATLAS 30 project, have similar morphological classifi¬ 
cations, have nearly identical surface brightness profiles 
(Sersic n ~ 3-4) and GC surface density profile shapes 
(Sersic n ~ 1.2-1.8). Yet, their GC systems show dif¬ 
ferent kinematics: VCC 1231’s GCs show little or no 
rotation. VCC 2000 GCs have a rotation amplitude 
consistent with the stars, but with a different rotation 
axis (PAkin)- VCC 1062’s GCs show fast rotation con¬ 
sistent with the stars both in amplitude and position 
angle. VCC 685’s GC system exhibits counter-rotation. 
This diversity suggests that these galaxies, despite their 
otherwise similar characteristics, had different assembly 
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Figure 8. Same as in Figure [5] but for VCC 2000 
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Figure 9. The kinematics fitting result as a function of radius in VCC 2000. The symbols are shown in the upper left. The SAURON 
data and GDL show rotation around the same axis, but the GDL data (which is at larger radius) has a smaller rotation amplitude than 
the SAURON data. The velocity dispersion profile of the GCs decreases with radius. The PAkin of the blue GCs is largely not coincident 
with the PAkin of the GDL. 


histories. Below, we discuss each of the galaxies in more 
detail. 

• VCC 1231 : Although the GC system as a whole 
does not show significant rotation, the red and blue 
GC individually appear to have rotation around 
nearly opposite kinematic axes, and with the red 
GCs aligned with the rotation of the stars. This 
finding is currently of low significance due to the 
low number of GCs in each sample. Many obser¬ 
vations suggest that the red GCs are more associ- 


ated with the stellar main body (Peng et aI1l2004at 
Faifer et all 120111 : iStrader et al. 120111 : iPota et al l 
20131) . It will be interesting to observe more GCs 
to test if the metal-poor GCs are truly counter¬ 
rotating with respect to the metal-rich GCs and 
the stars. If true, and whether the metal-poor 
GCs are accreted from low-mass progenitors (e.g., 
iCote et ahlll998D . then this galaxy requires a for¬ 
mation scenario where the angular momenta of the 
dry merged components and the gaseous compo¬ 
nents are nearly opposite. 
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VCC 2000 



Figure 10. Radial velocities versus position angle in VCC 2000. 
The GCs and the GDL both show rotation, but with different ro¬ 
tation axes. 

• VCC 2000 : This galaxy’s GC system is a more 
significant example of where the kinematics of blue 
GCs are different from those of the stars (see Figure 
HUD- In this case the rotation axis of the blue GCs 
is orthogonal to that of the stars. This galaxy GC 
system is consistent with a scenario where blue GCs 
are accreted from low-mass progenitor, and where 
the kinematics of these outer components are de¬ 
coupled from the material that eventually settles 
to form the stellar disk. 

• VCC 1062: This galaxy’s kinematics are the “sim¬ 
plest” in our sample. The red and blue GCs both 
have kinematics consistent with the stars. All three 
components of this galaxy show fast rotation about 
a common axis, the minor axis. In this case, the 
angular momenta of the progenitors that formed 
these components were all aligned. 

• VCC 685: The kinematics of the GCs in 
VCC 685 are the most complex and tantaliz¬ 
ing. Kinematically decoupled cores, or kinemat¬ 
ically distinct components (KDCs), have been 
found in the stellar kinematics in the central 
regions of many galaxies. This designation 
refers to an abrupt kinematic change between 


then the inner and outer re 
(IR.ix et al.! 1992 IMehlert et all 

gions ot a galaxy 
19981 IWernli et al. 

2002 Krainovic et all 20111 Arnold et alJ 2014 

Toloba et al. 201411. In VCC 

685, we find evi- 


dence that the GCs in the inner and outer regions 
are counter-rotating with respect to one another, 
with the inner component roughly aligned with the 
stars. This suggests discrete merging events. In¬ 
deed, for stellar KDCs, gas-rich mergers are often 
invoked as an explanation. In this case, however, 
the kinematic transition happens at much larger 
radii than for the stella r KDCs that are tr adition- 
ally discussed (although I Arnold et al.11'20 l~il find ev¬ 
idence of this behavior in some of their sample 
galaxies). More GC data, more stellar kinematic 
data at larger radii, and realistic simulations, are 
needed to investigate the counter-rotation in this 
galaxy. 


Lastly, we also analyze the kinematics of GCs subdi¬ 
vided by GC luminosity. Table [5] shows the rotation and 
velocity dispersion fitting results for subsamples of bright 
and faint GCs. We choose to divide the GC sample at 
the magnitude at which we have equal numbers of GCs 
in the two subsamples. The divided magnitude is shown 
in the last column of Table [3] For VCC 1231, we find 
that the brighter GCs have significantly more rotation 
than the fainter ones, which are consistent with no rota¬ 
tion. This is possibly because the brightest GCs are more 
likely to be blue, and hence their kinematics mirror the 
blue GCs. The fainter sample contains both blue and red 
GCs, and so their kinematics are a combination of these 
apparently counter-rotating populations. The other two 
galaxies do not show very significant differences in the ro¬ 
tation properties between the bright and faint GCs. VCC 
1062’s bright GCs, however, have a noticeably lower ve¬ 
locity dispersion (77± 17 km/s) compared to their fainter 
counterparts (134 ± 28 km/s). These difference worth 
further exploration in the future with larger data sets. 

5.2. Specific Angular Momentum 

To provide a more quantitative look at the kinematics 
of early-type galaxies, we use the specific angular mo¬ 
mentum. An early study of t he M87 GC system used a 
dime nsionless spin parameter (jKissler-Patig fc Gcbhardl] 
1998fl. We use t he di mensionless parameter defined by 
Emsellcm et all (I2007T ) to distinguish fast and slow rota¬ 
tors based on SAURON data within lR e . The expression 
for A is given by: 

A = (8) 

where Fi, Ri, Vi and a% are the flux, radius, velocity 
and velocity dispersion in the ith bin, respectively. Em- 
sellem et al. (2011) define the dividing line between slow 
and fast rotators at 

A = 0.31 x (9) 

where e is th e cllipticity of the system. However, 
lEmsellem et alJ (l2011l) ’s IFU study of ETC specific an¬ 
gular momentum was generally limited to stars within 
<lRe- 

Unlike for filled IFU observations, our GC and GDL 
data only sparsely sample the outer halos of our sam¬ 
ple galaxies. The density of our data, therefore, is not 
an exact representation of the true density of the GC 
or stars. Therefore, rather than sum over the full two- 
dimensional velocity field, which we do not have, we cal¬ 
culate A in one-dimension as follows: (1) Through our 
kinematic fitting, we obtain the rotation velocity and 
velocity dispersion profile with radius. This is done in 
elliptical annuli whose ellipticities are alternately set to 
two values: zero and the mean cllipticity of the stars. (2) 
For every bin in radius, we replace Fi in Equation |8] with 
the number of GCs expected at that radius, and use the 
rotation velocity at that radius. We call the A calculated 
in this way Aid. There is certainly a bias in using a 
one-dimensional method compared to A calculated with 
full two-dimensional data (which we now call A 2 d)- This 
bias is due to the fact that the 1-D method assigns the 
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Figure 11. Same as in Figure l5l but for VCC 1062. 
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Figure 12. The kinematic fitting results as a function of radius in VCC 1062. Red GCs, blue GCs, and the GDL all show significant 
rotation around the same rotation axis. The velocity dispersion profile of the red and blue GCs stays flat between 3 and 5 kpc, but the 
blue GC dispersion profile increases beyond 5 kpc. 
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VCC 1062 



Figure 13. Radial velocities versus position angle in VCC 1062. 
The GCs and GDL both show significant rotation about the same 
axis. 

peak velocity amplitude to every position at a given ra¬ 
dius. In reality, the velocity along the minor axis will be 
smaller than along the major axis. We can calibrate this 
bias, however, by comparing Aid and A 2 d for galaxies 
where we have full IFU coverage. To do this, we apply 
our one-dimensional method to the 260 ETG galaxies 
from ATLAS 30 . The result is shown in Figure [171 We 
can see that, as expected, Aid is generally overestimated 
compared to A 2 d- However, they follow a well-defined 
relation described as: 

A 2 d ~ 0.82 x Aid (10) 

A similar relation was also found between the 
SAURON IFU data an d long-slit ( V/a ) measurements 
(jOannellari et all 120071 ). and in t he long-slit studie s 
of A in early-type dwarf galaxies dToloba et al.1 120151) . 
Our scalin g factor is larger than the value found by 
iToloba et all (120151) due to the fact that all of their long- 
slit data is along the major axis. 

We apply this one-dimensional method to the GC data 
in our Gemini sample. We also include 6 galaxies with 
GC data from the literature (Table [6]). The uncertainties 
are determined from Monte Carlo simulations. Because 
the ellipticity of the GC system is not easily measured, 
we consider two cases, one where the ellipticity of the 
GC system is zero, and one where it is equal to that of 
the stars. The GC results are listed in Table [6] and Ta¬ 
bic [7] In practice, the difference is small between these 
two cases, and our results are not dependent on the as¬ 
sumed ellipticity of the GC system. In Figure HH we 
compare Aid of GC (blue circles) and SAURON data 
(red circles) with ellipticity. In all cases, we use the ellip¬ 
ticity of the stars. Although this is potentially an issue, 
studies of the GC systems of ETGs in this mass range 
show that there is some correlation between the overall 
shape of GC systems au d the ellipticity of the stellar dis¬ 
tribution (jWang et al.ll2013 ). The dashed lines link the 
same galaxy. The criterion of dividing slow and fast ro¬ 
tator with Equation [9] is rescaled by a factor 0.82 and 
plotted with a solid black line. 

The three stellar slow rotators (red points below the 
line), are massive ETGs (NGC 4486, 4472 and 4636), 
and their GCs also are relatively slowly rotating, with 


Xid,gc ~ 0-2- Massive ETGs may experience multiple 
mergers with different orbital configurations and mass ra¬ 
tio, and this leads to aggregate slow rotation in the halo. 
However, there are also cases of massive ETGs that have 
significant rotation in their outer regions even though the 
galaxies are classified as slow rotators (e .g. NGC 1407, 
Romanowskv et aO2009HPota et al.ll2013l and NGC 5128 
Peng et aLll2004bD . This case can be explained by angu¬ 
lar momentum transfer to the outer regions by mergers. 

For stellar fast rotators, however, the range in Aid.gc 
is large. The most obvious examples are our three inter¬ 
mediate luminosity galaxies, where Aid.gc of VCC 1062, 
2000 and 1231 is 0.9, 0.3 and 0.1, respectively. Half of 
stellar fast rotators show the Aid of GC are smaller than 
the Aid of star. Two stellar fast rotators (NGC 3379 
and 4473) have Aid.gc < 0-2, which means their outer 
regions are rotating very slowly. Similar results are also 
found by studying V ro t/cr (see the Figures 19 and 20 in 
iPota et ahll2013l) . It seems that, at least for the less mas¬ 
sive galaxies, GCs may be kinematically decoupled from 
the stars. 

5.3. Stellar kinematic trends with galactocentric radius 

ETGs at z ~ 2 appear to be smaller and more dense 
than present-day ETGs of similar mass dDaddi et ahl 
200,51 [Truiillo et al l l2007t iCenarro fc Truiillol 120091 
van Dokkum et al.l 120101 iWhitaker et al.l 120121) . A 
two phase mode l of galaxy assembly dNaab et al.1 120091 
lOser et all 120101 iKhochfar et al.1120111) has been used 
to explain these observations. In the “in situ” phase, 
dissipative processes such as gas-ri ch major mergers 
and fragmenting tu rbulent disks (iDekcl ct al.1 [2009i : 
iCeverino et al.ll2010D form the central region. In the “ex 
situ” phase, accreted stars from minor mergers expand 
a galaxy’s outskirts. These two phases may create dif¬ 
ferent kinematic components in ETGs, although it is not 
necessarily clear which component would have larger an¬ 
gular momentum. For example, the core region is likely 
to dissipate its angular momentum through violent relax¬ 
ation, but a late binar y merger can als o imprint rotation 
in the inner region s dCox et al.1120061 iBois et all 120111 
iMoodv et al.ll20ll) . The steady accretion processes in 
the outer region could create a dynamically hotter with 
low rotation, but could also imprint rotation if the ac¬ 
cretion is anisotropic. However, if there are indeed two 
regimes of galaxy building, we may expect to see this 
expressed in the kinematics. Indeed, some observations 
have found different stellar kinematics bet ween the inner 
and outer regions of early-type galaxies dProctor et al.l 
120091 iGoccato et aDl2009llArnold et al.ll20ll b 

One of the more comprehensive and rec ent studies of 
ETG stellar kinematics at large radii is by lArnold et al.l 
(|2014l) . They studied stellar kinematics in 22 nearby 
ETGs out to 2-4i? e with multi-s li t obs ervations from 
the SLUGGS survey dBrodie et al.l 120141) . They found 
that galaxies classified as fast rotators based on data in 
their inner regions can show a range of A profiles in their 
outer regions, suggesting that “the apparent unification 
of the elliptical and lenticular (SO) galaxy families in the 
ATLAS 30 survey may be a consequence of a limited field 
of view”. We measure Xid,gdl in the range of ~ 1.5- 
3 R e , and compare this measurement to Xm, sauron , 
which represents the specific ang ula r mom entum in the 
inner regions. As in lArnold et al.l (I2014D . we also find 
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Figure 15. The kinematics fitting result as a function of radius in VCC 685. The blue GCs exhibit counter-rotation, showing rotation in 
both the inner and outer regions, but around different rotation axes. The velocity dispersion profile of GC stays flat with radius. 


that A in the outer regions can differ from what is found 
in the inner regions. 

When trying to put into context the angular momen¬ 
tum profiles of galaxies, however, it is also important 
to consider ellipticity profiles. A changing ellipticity 
with radius leads to different expectations for a galaxy’s 
A with radius. For example, it is not s urprisin g that 
NGC 3377 shows a decrease in its A profile (jArnold et al.l 
I20l4h . as its stellar light also gets rounder with radius. 
In Figure 1191 we plot our measurements of Aid for the 
SAURON and GDL data versus radius, and also show the 
corresponding ellipticity profile (e), derived from NGVS 


imaging dFerrarese et al.1 [2012 1. For each galaxy’s Aid 
profile, the solid and dotted lines represent the SAURON 
and GDL data, respectively. Both lines are cumulative 
Aid profiles with radius, but each data set is treated inde¬ 
pendently to allow for comparison (i.e., the GDL profile 
does not include the SAURON data). 

Each of the three galaxies exhibits different character¬ 
istics. In VCC 1231 (left, red), Aid clearly decreases even 
though the ellipticity increases. This is not what one ex¬ 
pects for an anisotropic rotator at constant inclination. 
In VCC 2000 (center, green), the rotation declines as the 
galaxy gets rounder. In VCC 1062 (right, blue), the el- 
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Table 6 

General properties of literature galaxies 


ID 

R e 

M k 

D 

^star 

Ngc 

^lD,GC,e=0 

^1 D ,GC,e=star 

Rotation 

Ref. 


(kpc) 

(mag) 

(Mpc) 



< 1 Re 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

NGC 4494 

4.3 

-24.2 

16.6 

0.13 

117 

0.46 ± 0.04 

0.48 ±0.04 

F 

F±ll 

NGC 3379 

2.0 

-23.7 

10.2 

0.15 

52 

0.10 ± 0.03 

0.13 ± 0.03 

F 

P±04; P±06; B±06 

NGC 4636 

6.1 

-24.3 

14.2 

0.16 

459 

0.21 ± 0.02 

0.22 ± 0.02 

S 

S±12 

NGC 4472 (VCC 1226, M49) 

7.7 

-25.6 

16.7 

0.19 

263 

0.17 ± 0.04 

0.22 ± 0.05 

S 

C±03 

NGC 4649 (VCC 1978, M60) 

7.7 

-25.3 

16.3 

0.19 

121 

0.53 ±0.02 

0.56 ± 0.02 

F 

H±08 

NGC 4486 (VCC 1316, M87) 

13.1 

-25.3 

16.5 

0.14 

1004 

0.20 ± 0.01 

0.19 ±0.01 

S 

S+ll; Z±15* 


Note. — The effective radii (2) are from Ferrarese ct ah ( 20061) an d lEmsellem et al.| J20111) . The AT-band absolute ma gnitude (3) and dis tance 
(4) are fro m Pota ct al. (201^), Tables 1 and 7, where (3) is derived from 2MASS apparent magnitude and (4) are from [Tonrv ct al. (200l]) and 
IMei et all (1200 7h . (5) is the ellipticity with 1 — (h/a)^, where the axis ratio are from 2MASS. (6) is the total number of GCs with measured 
radial velocities. (7) and (8) are the projected specific angular momentum with e = 0 and e = Cstar • (9) Whether a galaxy is classified as a fast 
or slow rotator in ATLAS 3D lEmsellem_&t^_alJ12011J) . (10 ) Ref erences: F+iTjFostei^etal .1 (201 11 1 ^_P-j-04jPuziaet^al .1 1200^1 (JrM-06p|PiercejjyaI1 
1200611 : I3-t-06j|Bergond_et al. H2006l l: S+12 (fSc^iberth et ani2012T ): C-t-03 (ICote e^^T2003il : H+08 (IHwang et aflT2008lT : S-t-11 IStrader e^^T2^p]f : 
Z-tT5IZhang et^dll2015ll. 

* New GC data from our MMT observations, described by I Zhang et all (1201511 . are combined with previously published data. 

lipticity is nearly constant, but the rotation is increasing 
with radius. It is clear that while these galaxies may 
have been selected to be similar, their kinematics are 
quite different. 

Whether this supports a picture of “two-phase” for¬ 
mation is less clear. One expects continuous transitions 
from the inner to outer regions of galaxies that may 
not neces sarily result from disti nct mergers or formation 
episodes. iR.askutti et all (120141) have studied the stellar 
kinematics of 33 massive ETGs with IFU data out to 
2-5 R e and do not find a distinct transition between the 
inner and outer regions. Most of their sample galaxies 
are very massive, however, and it is possible that differ¬ 
ent kinematics between inner and outer regions are more 
easily seen in less massive systems, which are likely to 
have undergone fewer mergers than their massive coun¬ 
terparts. Individual merger events in massive galaxies 
are therefore more likely to be smoothed out by the larger 
number of merger and accretion events. 

6. CONCLUSIONS 

We have presented a kinematical analysis of halo stars 
and GCs in four intermediate luminosity galaxies in 
the Virgo cluster, based on Gemini/GMOS observations. 
These observations push these kinds of studies to lower 
mass ranges than previously done. Our galaxy diffuse 
light (GDL) and GC data extend out to 4R e and 8- 
14R e , respectively. We compare the GCs with the stars, 
and also compare the stellar kinematics in the inner 
(R < 1-R e ) and outer regions. Although the four galaxies 
in our sample are all uniformly classified as “fast rota¬ 
tors” by the ATLAS 3D project, we find that these galax¬ 
ies show a diversity of GC and stellar kinematics in their 
outer regions. 

1. The GCs display a wide range of kinematic prop¬ 
erties that in some galaxies follow the stars of their 
host galaxy, but in others do not. We see cases of 
nearly zero net rotation in the GCs (VCC 1231), 
rotation about a misaligned axis (VCC 2000), 
strong rotation (VCC 1062), and counter-rotation 
(VCC 685). 

2. We calculate the specific angular momentum, A, of 
the GCs for three galaxies in our sample as well as 


VCC 685 



Figure 16. Radial velocities versus position angle in VCC 685. 
The purple and green circles represent the inner and outer region 
GCs, respectively. The fits show strong counter-rotation. The ma¬ 
genta curve and points represents the rotation of the GDL, which 
lines up more with the rotation of the inner GCs. 



Figure 17. The comparison of Ai£> and A 2 d- The black points 
and the red points represent the slow and fast rotator. The black 
line is a 1:1 relation. The blue line show the best fitting linear 
relation to the data. 
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Figure 18. The specific angular momentum of GCs and SAURON data as function of ellipticity . The red circles represent the X\d of 
the GC. The blue circles represent the X\d of the SAURON data. The dash lines link the same galaxy. The solid line is Xid = 0.38 x y/e. 
Most stellar slow rotators also show slow rotation in their GC systems, with Xqc ~ 0.2. Stellar fast rotators display a large range in A gc- 


Table 7 

Specific angular momentum 


ID 

(1) 

Vo 

(SAURON) 

(2) 

Al D 

(SAURON) 

(3) 

Aid 

(GDL) 

(4) 

AlD,GC,e=0 

(GCs) 

(5) 

Al D,GC,e=star 

(GCs) 

(6) 

VCC1231 

0.23 

0.34 

0.17±0.02 

0.04±0.04 

0.08 ± 0.06 

VCC2000 

0.55 

0.70 

0.49±0.02 

0.30±0.06 

0.32 ±0.06 

VCC1062 

0.36 

0.51 

0.84±0.01 

0.82±0.03 

0.88 ±0.02 


Note. — (2) A 2 d measurement from SAURON data, (3) \\d measurement from 
SAURON data, (4) X\d measurement from Gemini (GDL), (5) and (6) X\d,gc f° r the 
GCs assuming different ellipticities. The uncertainties for the SAURON data are very 
small, and they are not listed in the table. 

in 6 galaxies from the literature. We find that mas¬ 
sive ETGs that are classified as stellar slow rotators 
in their inner R < 1 R e by ATLAS 30 also tend to 
have slowly rotating GC systems (although at least 
one counter-example exists in the literature). The 
GC systems of stellar fast rotators, however, dis¬ 
play a wide range of A. Some galaxies whose stars 
are rapidly rotating have slowly or only moderately 
rotating GC systems, and vice versa. We interpret 
this diversity as the consequence of the stochastic 
nature of merging processes in less massive galax¬ 
ies. 

3. We compare the specific angular momentum for 
stars in the inner (R < 1 R e ) and outer (1.5 < 

R < 3.0) regions of three of our sample galaxies. 


In conjunction with the ellipticity profiles derived 
from deep imaging, we find that one of our galax¬ 
ies (VCC 1231) has a A profile that decreases de¬ 
spite increasing ellipticity with radius, which is not 
what is expected for a simple anisotropic rotator 
at fixed inclination. The other two galaxies show 
falling (VCC 2000) and rising (VCC 1062) A pro¬ 
files. Even in this small sample, we can see that 
the stellar kinematics of otherwise morphologically 
similar galaxies can exhibit a diversity of properties 
at larger radii. 

The kinematic differences we and others have found 
between GCs and stars, and between the inner and outer 
regions, in intermediate luminosity ETGs are beginning 
to allow us to dissect the merging and accretion histories 
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Figure 19. The Ai £> and the cllipticity profiles as radius for VCC 1231, 2000 and 1062. In the top panel, the solid curve represents the 
Ai D profile measurements from SAURON data. The dash curve represents the Aio profile measurement from the GDL. The individual 
A | D profiles (SAURON and GDL) are each cumulative with radius, but the two data sets are treated separately to allow for a comparison. 


of normal galaxies. The complex array of properties we 
are finding shows that this kind of wide-held kinematic 
study of a larger sample of intermediate-mass galaxies 
will be needed to fully sample the parameter space of 
galaxy assembly. 
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Figure 20. The 2D velocity and velocity dispersion fiedls of VCC 1231 and 2000, without smooth and double data. 


APPENDIX 

UNSMOOTHED AND RAW VELOCITY FIELD 

The unsmoothed and raw GC and stellar (GDL) velocity and velocity dispersion fields are shown in Figure [20] and 
Figure [21]. The rotation information we get from these figures is consistent with the smoothed and point symmetric 
velocity fields in Figures [5] ® HU and [14] . Our quantitative analysis of specific angular momentum is based on the 
unsmoothed and non-symmetrized data. 

The 2D velocity dispersio n fiel d of VCC 123 1 sh ows good agreement with SAURON in the overlap region. Our 
result is also consistent with lFoster et all (12013T) and lArnold et al.1 (I2014T) . who find two velocity dispersion peaks at 3 
kpc along the major axis. VCC 1231 is also called as “2cr” galaxy. 
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Figure 21. 
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Same as in Figure [201 but for VCC 1062 and 685. 























Table 8 

Globular Cluster 
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ID 

R.A.(J2000) 

Decl.(J2000) 

^OS 

9 

z 

9-z 

E(B-V) 


(deg) 

(deg) 

(km s _1 ) 

(mag) 

(mag) 

(mag) 

(mag) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

VCC 1231 

VCC1231GC-01 

187.420170 

13.436530 

24504=19 

21.2804=0.016 

20.229=10.024 

0.9914=0.007 

0.028 

VCC1231GC-02 

187.422675 

13.437170 

2107d=22 

21.9604=0.026 

20.9224=0.029 

0.9734=0.033 

0.028 

VCC1231GC-03 

187.428810 

13.436997 

23774=20 

22.1924=0.022 

20.8934=0.020 

1.2534=0.028 

0.028 

VCC1231GC-04 

187.413315 

13.420913 

2429±15 

22.2334=0.030 

20.7884=0.019 

1.3744=0.024 

0.028 

VCC1231GC-05 

187.418970 

13.424147 

23574=17 

23.0334=0.032 

22.1964=0.031 

0.8174=0.022 

0.028 

VCC1231GC-06 

187.411215 

13.419499 

22294=21 

23.1474=0.038 

22.1854=0.031 

0.9074=0.025 

0.028 

VCC1231GC-07 

187.405065 

13.453014 

20074=18 

22.0964=0.007 

21.3594=0.016 

0.7034=0.012 

0.028 

VCC1231GC-08 

187.429815 

13.401846 

20584=49 

22.9334=0.013 

22.2714=0.037 

0.5514=0.030 

0.028 

VCC1231GC-09 

187.432710 

13.444921 

22824=20 

22.8554=0.022 

21.7904=0.021 

1.0184=0.025 

0.028 

VCC1231GC-10 

187.441245 

13.429025 

24354=28 

21.8424=0.037 

20.9774=0.018 

0.7994=0.039 

0.028 

VCC1231GC-11 

187.437780 

13.422738 

25994=24 

21.8214=0.012 

20.9034=0.020 

0.8324=0.019 

0.028 

VCC1231GC-12 

187.447410 

13.432945 

23244=20 

21.2094=0.029 

20.1544=0.020 

0.9924=0.028 

0.028 

VCC1231GC-13 

187.450110 

13.433290 

25334=21 

21.8134=0.019 

20.8174=0.012 

0.9464=0.020 

0.028 

VCC1231GC-14 

187.454700 

13.403273 

22504=24 

22.4854=0.018 

21.5704=0.016 

0.8714=0.022 

0.028 

VCC1231GC-15 

187.452990 

13.445077 

23564=20 

21.8614=0.017 

20.3434=0.013 

1.4364=0.016 

0.028 

VCC1231GC-16 

187.456470 

13.424834 

23894=27 

22.0054=0.020 

20.9624=0.015 

0.9904=0.021 

0.028 

VCC1231GC-17 

187.463295 

13.427741 

23674=17 

21.3064=0.014 

20.3204=0.017 

0.9204=0.019 

0.028 

VCC1231GC-18 

187.474695 

13.439762 

21514=26 

22.3594=0.022 

21.4124=0.016 

0.8854=0.022 

0.028 

VCC1231GC-19 

187.468155 

13.427391 

25234=17 

21.8064=0.032 

20.4694=0.011 

1.2744=0.028 

0.028 

VCC1231GC-20 

187.475715 

13.426306 

21074=29 

22.5504=0.020 

21.5214=0.019 

0.9674=0.022 

0.028 

VCC1231GC-21 

187.469025 

13.435520 

23374=20 

22.9324=0.025 

21.6074=0.014 

1.2784=0.024 

0.028 

VCC1231GC-22 

187.479150 

13.444892 

21484=23 

22.6734=0.023 

21.6394=0.015 

1.0144=0.026 

0.028 

VCC1231GC-23 

187.485810 

13.446020 

24124=15 

21.1664=0.032 

20.1944=0.022 

0.8804=0.022 

0.028 

VCC1231GC-24 

187.496295 

13.436360 

22004=11 

21.1834=0.004 

20.4674=0.008 

0.6474=0.007 

0.028 

VCC1231GC-25 

187.429695 

13.433433 

19094=15 

22.6674=0.017 

21.2184=0.019 

1.3844=0.019 

0.028 

VCC1231GC-26 

187.424715 

13.427611 

19884=23 

22.6164=0.017 

21.2614=0.018 

1.2944=0.019 

0.028 

VCC1231GC-27 

187.422900 

13.414083 

22264=26 

22.7854=0.056 

21.7374=0.062 

0.9094=0.021 

0.028 

VCC1231GC-28 

187.422360 

13.402969 

23644=14 

21.2744=0.003 

20.4454=0.007 

0.7634=0.006 

0.028 

VCC1231GC-29 

187.450515 

13.431453 

21984=31 

21.2064=0.014 

20.2874=0.018 

0.8494=0.014 

0.028 

VCC1231GC-30 

187.447620 

13.429241 

24184=23 

21.8434=0.025 

20.5144=0.016 

1.2624=0.027 

0.028 

VCC1231GC-31 

187.442595 

13.426672 

21514=20 

21.4114=0.016 

20.4674=0.013 

0.8764=0.015 

0.028 

VCC1231GC-32 

187.461765 

13.411035 

23514=16 

22.1394=0.080 

21.0534=0.022 

0.9994=0.074 

0.028 

VCC1231GC-33 

187.434585 

13.434356 

23974=16 

22.9534=0.027 

21.8844=0.025 

0.9904=0.025 

0.028 

VCC1231GC-34 

187.444185 

13.431620 

21554=69 

22.6714=0.023 

21.7584=0.023 

0.8264=0.023 

0.028 

VCC1231GC-35 

187.457610 

13.433330 

20134=24 

22.2744=0.020 

21.2764=0.020 

0.9314=0.023 

0.028 

VCC1231GC-36 

187.452285 

13.459685 

20124=38 

22.9524=0.034 

21.8524=0.026 

1.0314=0.031 

0.028 

VCC1231GC-37 

187.460490 

13.472438 

21404=21 

21.3614=0.005 

20.3004=0.007 

1.0464=0.006 

0.028 

VCC1231GC-38 

187.473540 

13.428547 

20334=35 

22.9644=0.053 

22.0824=0.042 

0.8214=0.064 

0.028 

VCC1231GC-39 

187.493505 

13.426348 

19294=34 

20.8094=0.002 

19.9554=0.005 

0.7554=0.005 

0.028 

VCC1231GC-40 

187.428210 

13.420699 

23714=32 

23.1314=0.017 

21.8014=0.012 

1.2714=0.019 

0.028 

VCC1231GC-41 

187.429515 

13.423301 

21994=48 

23.3474=0.025 

22.4044=0.022 

0.8854=0.029 

0.028 

VCC1231GC-42 

187.451670 

13.442974 

25384=16 

22.3474=0.016 

21.1814=0.025 

1.1144=0.029 

0.028 

VCC1231GC-43 

187.457820 

13.406286 

22724=15 

21.5644=0.016 

20.3854=0.018 

1.0754=0.014 

0.028 

VCC1231GC-44 

187.465260 

13.417413 

24674=11 

20.9554=0.024 

19.7764=0.012 

1.1264=0.023 

0.028 

VCC1231GC-45 

187.453365 

13.425970 

22304=14 

21.3044=0.021 

20.2724=0.015 

0.9574=0.017 

0.028 

VCC1231GC-46 

187.438560 

13.435633 

20144=17 

22.7614=0.015 

21.6144=0.020 

1.0704=0.020 

0.028 

VCC1231GC-47 

187.445190 

13.432550 

22134=36 

23.1514=0.026 

21.7274=0.022 

1.3584=0.031 

0.028 

VCC1231GC-48 

187.461540 

13.434665 

21254=27 

22.7844=0.024 

21.8864=0.016 

0.8374=0.026 

0.028 

VCC1231GC-49 

187.474830 

13.443451 

24504=20 

22.7304=0.021 

21.4254=0.020 

1.2364=0.026 

0.028 

VCC1231GC-50 

187.468185 

13.430573 

24594=38 

23.1874=0.015 

21.9054=0.026 

1.2064=0.021 

0.028 

VCC1231GC-51 

187.476593 

13.420838 

21614=24 

22.7404=0.012 

21.8864=0.028 

0.6784=0.021 

0.028 


VCC 2000 


VCC2000GC-01 

191.092740 

11.213082 

10814=12 

20.5724=0.002 

19.8934=0.005 

0.6114=0.004 

0.034 

VCC2000GC-02 

191.122215 

11.219080 

9494=42 

22.7414=0.027 

21.9124=0.027 

0.7784=0.025 

0.034 

VCC2000GC-03 

191.125230 

11.215344 

9074=28 

22.4784=0.029 

21.5564=0.029 

0.9064=0.022 

0.034 

VCC2000GC-04 

191.128680 

11.213797 

10164=19 

22.7274=0.066 

21.8024=0.052 

0.8894=0.028 

0.034 

VCC2000GC-05 

191.109285 

11.205482 

9204=17 

21.1514=0.015 

20.2434=0.014 

0.8894=0.022 

0.034 

VCC2000GC-06 

191.138010 

11.200358 

10744=24 

22.5484=0.069 

21.2584=0.023 

1.1814=0.070 

0.034 

VCC2000GC-07 

191.135130 

11.198518 

7894=22 

21.8854=0.032 

20.8824=0.063 

0.8634=0.032 

0.034 

VCC2000GC-08 

191.131050 

11.197987 

9174=33 

22.5234=0.032 

21.5464=0.013 

0.8954=0.035 

0.034 

VCC2000GC-09 

191.114865 

11.196526 

10094=18 

22.7784=0.027 

21.6424=0.014 

1.0714=0.028 

0.034 

VCC2000GC-10 

191.149545 

11.201708 

9904=25 

22.2504=0.023 

21.3054=0.020 

0.8884=0.025 

0.034 

VCC2000GC-11 

191.143635 

11.196664 

13444=48 

22.4564=0.036 

21.6224=0.029 

0.7604=0.031 

0.034 

VCC2000GC-12 

191.140995 

11.196223 

12634=42 

22.2474=0.019 

21.4274=0.014 

0.7524=0.020 

0.034 

VCC2000GC-13 

191.146065 

11.194601 

11334=17 

22.1564=0.019 

20.9934=0.015 

1.1164=0.021 

0.034 

VCC2000GC-14 

191.151405 

11.193774 

9724=78 

22.5614=0.020 

21.8074=0.013 

0.7054=0.022 

0.034 

VCC2000GC-15 

191.167200 

11.189954 

10854=28 

21.3894=0.004 

20.7584=0.011 

0.5764=0.007 

0.034 

VCC2000GC-16 

191.119545 

11.199561 

11714=22 

22.4624=0.017 

21.4554=0.015 

0.9594=0.019 

0.034 

VCC2000GC-17 

191.125890 

11.190406 

12504=18 

22.2114=0.033 

20.8464=0.015 

1.2794=0.025 

0.034 

VCC2000GC-18 

191.107575 

11.187544 

9824=31 

22.3564=0.015 

21.5314=0.021 

0.7424=0.021 

0.034 

VCC2000GC-19 

191.123100 

11.191837 

12674=15 

21.0294=0.024 

20.0324=0.022 

0.9074=0.021 

0.034 
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Table 8 — Continued 


ID 

(1) 

R.A.(J2000) 

(deg) 

(2) 

Decl.(J2000) 

(deg) 

(3) 

^los 

(km s -1 ) 

(4) 

9 

(mag) 

(5) 

z 

(mag) 

(6) 

9-z 

(mag) 

(7) 

E(B-V) 

(mag) 

(8) 

VCC2000GC-20 

191.129370 

11.193860 

12704=19 

21.6114=0.027 

20.6324=0.019 

0.9034=0.011 

0.034 

VCC2000GC-21 

191.111385 

11.195764 

12144=43 

23.0324=0.029 

22.2344=0.030 

0.7194=0.034 

0.034 

VCC2000GC-22 

191.141985 

11.195713 

9214=27 

22.4674=0.032 

21.4664=0.018 

0.9424=0.033 

0.034 

VCC2000GC-23 

191.138505 

11.193366 

841±17 

21.2664=0.019 

20.2634=0.013 

0.9334=0.017 

0.034 

VCC2000GC-24 

191.146710 

11.191358 

10224=18 

20.5534=0.058 

19.6914=0.052 

0.7874=0.023 

0.034 

VCC2000GC-25 

191.155815 

11.187847 

9824=23 

21.8484=0.029 

20.9554=0.023 

0.8704=0.031 

0.034 

VCC2000GC-26 

191.150445 

11.183163 

9424=21 

21.5434=0.031 

20.7434=0.034 

0.7314=0.011 

0.034 

VCC2000GC-27 

191.145060 

11.197882 

8634=28 

22.8774=0.024 

22.0854=0.020 

0.7294=0.026 

0.034 

VCC2000GC-28 

191.159205 

11.174632 

11544=24 

21.1524=0.003 

20.4394=0.008 

0.6354=0.006 

0.034 

VCC2000GC-29 

191.138430 

11.210131 

7604=32 

22.5794=0.030 

21.7434=0.021 

0.7684=0.026 

0.034 

VCC2000GC-30 

191.127015 

11.186158 

10954=29 

22.5104=0.032 

21.3584=0.017 

1.0794=0.024 

0.034 

VCC2000GC-31 

191.119665 

11.189951 

8454=16 

20.7774=0.022 

19.7714=0.019 

0.9374=0.022 

0.034 

VCC2000GC-32 

191.130135 

11.184649 

1261±19 

21.5724=0.020 

20.5674=0.008 

0.9004=0.019 

0.034 

VCC2000GC-33 

191.133315 

11.183667 

1010=t27 

22.5944=0.032 

21.5634=0.024 

0.9704=0.036 

0.034 

VCC2000GC-34 

191.110065 

11.210609 

11124=33 

22.8394=0.022 

22.0134=0.020 

0.8254=0.023 

0.034 

VCC2000GC-35 

191.122260 

11.178163 

10874=28 

22.9254=0.038 

21.7794=0.016 

1.0834=0.035 

0.034 

VCC2000GC-36 

191.116050 

11.230333 

11094=20 

20.5994=0.002 

19.8784=0.006 

0.6294=0.004 

0.034 

VCC2000GC-37 

191.140260 

11.163837 

11894=26 

22.3804=0.012 

21.4394=0.014 

0.9024=0.016 

0.034 

VCC2000GC-38 

191.138190 

11.176321 

10044=38 

21.1934=0.030 

20.3224=0.026 

0.7924=0.017 

0.034 

VCC2000GC-39 

191.147055 

11.198483 

11494=43 

23.0004=0.026 

22.1694=0.025 

0.7774=0.023 

0.034 

VCC2000GC-40 

191.146830 

11.172589 

9704=56 

23.0304=0.033 

22.1744=0.033 

0.8094=0.024 

0.034 

VCC2000GC-41 

191.155275 

11.201397 

10604=14 

21.3654=0.015 

20.4654=0.011 

0.8514=0.013 

0.034 

VCC2000GC-42 

191.161890 

11.205892 

10864=22 

22.7104=0.037 

21.8154=0.026 

0.8434=0.031 

0.034 

VCC2000GC-43 

191.167455 

11.181426 

9424=37 

22.5614=0.010 

21.8124=0.028 

0.6584=0.019 

0.034 


VCC 1062 


VCC1062GC-01 

187.040970 

9.794455 

683=1=12 

22.890=10.020 

21.5344=0.020 

1.3254=0.023 

0.022 

VCC1062GC-02 

187.000530 

9.817132 

444=1=19 

21.895=t0.013 

20.9644=0.013 

0.8864=0.014 

0.022 

VCC1062GC-03 

187.004025 

9.816408 

512=t23 

22.575=t0.045 

21.7574=0.062 

0.7294=0.060 

0.022 

VCC1062GC-04 

186.995685 

9.807220 

306=1=18 

24.395=t0.038 

23.1274=0.057 

1.2314=0.059 

0.022 

VCC1062GC-05 

186.999240 

9.805346 

518=1=15 

21.978=t0.024 

21.0204=0.018 

0.9024=0.022 

0.022 

VCC1062GC-06 

186.989805 

9.792907 

3424=21 

22.503=t0.029 

21.5684=0.029 

0.8594=0.023 

0.022 

VCC1062GC-07 

186.991335 

9.802511 

2294=30 

23.130=t0.020 

22.1344=0.022 

0.9804=0.023 

0.022 

VCC1062GC-08 

186.985650 

9.798775 

3264=21 

22.977=t0.016 

21.8234=0.026 

1.1094=0.020 

0.022 

VCC1062GC-09 

187.012260 

9.817571 

421:121 

22.418=t0.012 

21.5704=0.025 

0.7914=0.017 

0.022 

VCC1062GC-10 

187.009785 

9.810227 

4174=21 

22.582=t0.021 

21.1964=0.017 

1.3054=0.023 

0.022 

VCC1062GC-11 

187.017090 

9.809849 

619=tl8 

22.279=t0.023 

21.2954=0.017 

0.9084=0.020 

0.022 

VCC1062GC-12 

187.026240 

9.811627 

595=t23 

21.870=t0.032 

20.9284=0.051 

0.8054=0.021 

0.022 

VCC1062GC-13 

187.025115 

9.796717 

679=t20 

21.588=t0.025 

20.7164=0.038 

0.7784=0.017 

0.022 

VCC1062GC-14 

187.037430 

9.792029 

516=tl3 

21.699=t0.017 

20.7404=0.012 

0.9144=0.020 

0.022 

VCC1062GC-15 

187.030080 

9.811315 

547=tl8 

22.869=t0.026 

21.9444=0.024 

0.8914=0.021 

0.022 

VCC1062GC-16 

187.047585 

9.805256 

583=tl7 

22.437=t0.021 

21.3094=0.016 

1.0884=0.021 

0.022 

VCC1062GC-17 

187.044615 

9.816442 

3014=25 

22.926=t0.020 

21.9794=0.015 

0.9004=0.022 

0.022 

VCC1062GC-18 

187.054260 

9.809342 

645=t21 

23.1024=0.038 

21.6814=0.027 

1.4134=0.025 

0.022 

VCC1062GC-19 

187.000560 

9.796535 

380=tl9 

22.4584=0.039 

20.9654=0.037 

1.4874=0.020 

0.022 

VCC1062GC-20 

187.005885 

9.798541 

472=tl7 

22.2544=0.049 

20.6624=0.061 

1.4824=0.033 

0.022 

VCC1062GC-21 

186.994635 

9.807293 

309=t21 

22.5744=0.020 

21.1134=0.018 

1.4144=0.022 

0.022 

VCC1062GC-22 

187.003020 

9.802365 

279=tll 

21.2284=0.017 

20.0094=0.020 

1.1654=0.019 

0.022 

VCC1062GC-23 

186.999375 

9.790121 

234=tl8 

23.1344=0.033 

22.1054=0.017 

1.0384=0.030 

0.022 

VCC1062GC-24 

186.986835 

9.806428 

311=114 

22.9794=0.016 

21.8924=0.028 

1.0944=0.019 

0.022 

VCC1062GC-25 

187.014075 

9.808210 

461=tl9 

22.6654=0.033 

21.2344=0.024 

1.2934=0.017 

0.022 

VCC1062GC-26 

187.018665 

9.791675 

486=tl3 

21.2684=0.019 

20.3284=0.015 

0.8594=0.019 

0.022 

VCC1062GC-27 

187.012890 

9.792312 

310=tl8 

22.9464=0.019 

21.5734=0.021 

1.2894=0.020 

0.022 

VCC1062GC-28 

187.028310 

9.813976 

605=tl9 

22.5924=0.020 

21.5474=0.017 

0.9914=0.018 

0.022 

VCC1062GC-29 

187.037685 

9.798342 

743=tl9 

21.7714=0.058 

20.9114=0.051 

0.8664=0.016 

0.022 

VCC1062GC-30 

187.025085 

9.798357 

488=t57 

23.1134=0.030 

22.1454=0.051 

0.8664=0.037 

0.022 

VCC1062GC-31 

187.050615 

9.787901 

648=tl7 

22.4094=0.009 

21.6214=0.021 

0.7184=0.015 

0.022 

VCC1062GC-32 

187.067940 

9.802051 

567=tll 

21.5064=0.004 

20.7494=0.010 

0.7204=0.007 

0.022 

VCC1062GC-33 

187.004535 

9.801900 

281=112 

24.7774=0.062 

23.4604=0.066 

1.2714=0.068 

0.022 

VCC1062GC-34 

187.011405 

9.804032 

401=tll 

21.1894=0.016 

20.0914=0.020 

1.0124=0.016 

0.022 

VCC1062GC-35 

187.018740 

9.803278 

569=t24 

21.8044=0.017 

20.6024=0.019 

1.1574=0.021 

0.022 

VCC1062GC-36 

187.026120 

9.804276 

717=tl8 

22.3994=0.021 

21.3254=0.021 

1.0024=0.024 

0.022 

VCC1062GC-37 

187.023465 

9.805268 

414=t23 

22.9744=0.022 

21.4964=0.045 

1.3664=0.030 

0.022 

VCC1062GC-38 

187.036740 

9.806686 

757=tl5 

22.8424=0.027 

21.2574=0.025 

1.5344=0.033 

0.022 

VCC1062GC-39 

187.001925 

9.804249 

336=tl0 

20.5264=0.020 

19.3574=0.023 

1.1114=0.023 

0.022 

VCC1062GC-40 

186.992850 

9.808851 

287=tl8 

22.1794=0.008 

21.1394=0.014 

0.9914=0.011 

0.022 


VCC 685 


VCC685GC-01 

185.989590 

16.651251 

13414=24 

22.2624=0.008 

21.6164=0.026 

0.6094=0.017 

0.028 

VCC685GC-02 

185.959680 

16.679079 

13344=27 

22.1304=0.021 

21.2054=0.014 

0.8684=0.019 

0.028 

VCC685GC-03 

185.975325 

16.684580 

14224=22 

21.9544=0.014 

20.9984=0.011 

0.8974=0.016 

0.028 

VCC685GC-04 

185.978040 

16.679291 

12854=31 

22.2564=0.014 

21.4224=0.017 

0.7544=0.017 

0.028 

VCC685GC-05 

185.988135 

16.683411 

9184=42 

22.3674=0.021 

21.2384=0.022 

1.0214=0.016 

0.028 
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Table 8 — Continued 


ID 

(1) 

R.A.(J2000) 

(deg) 

(2) 

Decl.(J2000) 

(deg) 

(3) 

t^los 

(km s -1 ) 

(4) 

9 

(mag) 

(5) 

z 

(mag) 

(6) 

9-z 

(mag) 

(7) 

E(B-V) 

(mag) 

(8) 

VCC685GC-06 

185.985465 

16.692770 

1025=1=23 

21.566=1=0.019 

20.658=1=0.017 

0.844=1=0.010 

0.028 

VCC685GC-07 

185.989155 

16.688009 

1062=1=28 

21.987=1=0.026 

21.125=1=0.018 

0.811=1=0.024 

0.028 

VCC685GC-08 

185.989965 

16.701780 

1168±14 

21.4844=0.021 

20.311=1=0.013 

1.115=1=0.023 

0.028 

VCC685GC-09 

185.991735 

16.694599 

1588±17 

21.529±0.012 

20.532=1=0.018 

0.903=1=0.017 

0.028 

VCC685GC-10 

185.994150 

16.696199 

1520=1=38 

21.916=1=0.023 

20.950=1=0.028 

0.910=1=0.018 

0.028 

VCC685GC-11 

185.991225 

16.709190 

1343=1=21 

20.597=1=0.018 

19.737=1=0.012 

0.811=1=0.006 

0.028 

VCC685GC-12 

185.997015 

16.701620 

1178=1=22 

21.662=1=0.016 

20.724=1=0.012 

0.878=1=0.020 

0.028 

VCC685GC-13 

185.996025 

16.711010 

1208=1=15 

21.145=1=0.015 

19.992=1=0.012 

1.084=1=0.017 

0.028 

VCC685GC-14 

186.006900 

16.716101 

925=1=33 

22.649=1=0.019 

21.725=1=0.015 

0.871=1=0.021 

0.028 

VCC685GC-15 

186.000240 

16.730289 

1098=1=31 

21.078=1=0.003 

20.493=1=0.012 

0.508=1=0.008 

0.028 


Note. — Col.(l): Object ID; Col.(2): Right ascension in decimal degrees (J2000); Col.(3): Declination in decimal degrees (J2000); Col.(4): 
Heliocentric radial velocity and uncertainty; Cols.(5,6): g and z band magnitudes from the ACSVCS and NGVS (not corrected for Galactic 
extinction); Col.(7): g — z color; Col.(8): The Galactic reddening determined by Schlegel et al. (1998). The corrections for foreground reddening 
were taken to be A g = 3.634 E(B — V ) and A z = 1A85E(B — V ) in the g and z bands, respectively dJordan et al.ir2004h - 


Table 9 

Galaxy Diffuse Light 


ID 

R.A.(J2000) 

Decl.(J2000) 

^los 

°fos 


(deg) 

(deg) 

(km s -1 ) 

(km s -1 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

VCC 1231 

VCC1231GDL-01 

187.450515 

13.431453 

2188=1=14 

174=1= 4 

VCC1231GDL-02 

187.447620 

13.429241 

2205=tl8 

208=1= 4 

VCC1231GDL-03 

187.453365 

13.425970 

2244±12 

153=1= 4 

VCC1231GDL-04 

187.456470 

13.424834 

2270=1=14 

157=1= 5 

VCC1231GDL-05 

187.457610 

13.433330 

2252=1=14 

154=1= 5 

VCC1231GDL-06 

187.447410 

13.432945 

2193=1=14 

185=1= 5 

VCC1231GDL-07 

187.450110 

13.433290 

2191=1=14 

163=1= 3 

VCC1231GDL-08 

187.445190 

13.432550 

2206=tl4 

189=1= 2 

VCC1231GDL-09 

187.444185 

13.431620 

2200±18 

197=1= 5 

VCC1231GDL-10 

187.441245 

13.429025 

2215=1=18 

204=1= 5 

VCC1231GDL-11 

187.442595 

13.426672 

2209=1=16 

195=1= 5 

VCC1231GDL-12 

187.461540 

13.434665 

2251=1=16 

165=1= 6 

VCC1231GDL-13 

187.463295 

13.427741 

2250=1=15 

203=1= 3 

VCC1231GDL-14 

187.468185 

13.430573 

2244=1=17 

197=1= 6 

VCC1231GDL-15 

187.468155 

13.427391 

2249=1=17 

196=1= 5 

VCC1231GDL-16 

187.473540 

13.428547 

2225=1=25 

186=1=11 

VCC1231GDL-17 

187.438560 

13.435633 

2229=1=22 

211=1=10 

VCC1231GDL-18 

187.475715 

13.426306 

2243=1=29 


VCC1231GDL-19 

187.469025 

13.435520 

2239=1=25 


VCC1231GDL-20 

187.434585 

13.434356 

2203=1=17 


VCC1231GDL-21 

187.437780 

13.422738 

2217=1=24 


VCC1231GDL-22 

187.428810 

13.436997 

2257=1=29 


VCC1231GDL-23 

187.429695 

13.433433 

2200=1=20 


VCC1231GDL-24 

187.429515 

13.423301 

2206=1=24 


VCC1231GDL-25 

187.451670 

13.442974 

2210=1=26 


VCC1231GDL-26 

187.465260 

13.417413 

2265=1=32 


VCC 2000 

VCC2000GDL-01 

191.134440 

11.189226 

1032=1=12 

142=1= 3 

VCC2000GDL-02 

191.131935 

11.189683 

1144=1=12 

145=1= 3 

VCC2000GDL-03 

191.138505 

11.193366 

1034=1=10 

138=1= 7 

VCC2000GDL-04 

191.129370 

11.193860 

1115=1=10 

129=1= 6 

VCC2000GDL-05 

191.125890 

11.190406 

1166=tll 

138=1= 9 

VCC2000GDL-06 

191.123100 

11.191837 

1139=1=12 

149=1=12 

VCC2000GDL-07 

191.127015 

11.186158 

1109=1:13 

115=1=10 

VCC2000GDL-08 

191.130135 

11.184649 

1094=1=16 

126=1=10 

VCC2000GDL-09 

191.133315 

11.183667 

1065=1=16 

118=1= 8 

VCC2000GDL-10 

191.138505 

11.193366 

1034=1= 9 


VCC2000GDL-11 

191.119665 

11.189951 

1134=1=17 


VCC2000GDL-12 

191.146710 

11.191358 

996=1=28 


VCC2000GDL-13 

191.141985 

11.195713 

1067=1=21 


VCC2000GDL-14 

191.140995 

11.196223 

1058=1=20 


VCC2000GDL-15 

191.135130 

11.198518 

1057=1=16 


VCC2000GDL-16 

191.131050 

11.197987 

1075=1=19 


VCC2000GDL-17 

191.143635 

11.196664 

1040=1=25 


VCC2000GDL-18 

191.146065 

11.194601 

1029=1=26 


VCC2000GDL-19 

191.145060 

11.197882 

995=1=33 


VCC2000GDL-20 

191.119545 

11.199561 

1052=1=44 
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Table 9 — Continued 


ID 

R.A.(J2000) 

Decl.(J2000) 

'^los 

^los 


(deg) 

(deg) 

(km s -1 ) 

(km s -1 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

VCC2000GDL-21 

191.119845 

11.183841 

1080=1=42 


VCC 1062 

VCC1062GDL-01 

187.018740 

9.803278 

611=1=17 

172=t 3 

VCC1062GDL-02 

187.026120 

9.804276 

683=1=11 

128=1= 6 

VCC1062GDL-03 

187.023465 

9.805268 

640=1=12 

123=1= 4 

VCC1062GDL-04 

187.014075 

9.808210 

486=1=11 

110=1= 5 

VCC1062GDL-05 

187.011405 

9.804032 

419d=12 

136=1=20 

VCC1062GDL-06 

187.017090 

9.809849 

531±12 

116=1= 7 

VCC1062GDL-07 

187.009785 

9.810227 

468=1=14 

115=1=10 

VCC1062GDL-08 

187.026240 

9.811627 

613±19 

133=1=12 

VCC1062GDL-09 

187.030080 

9.811315 

640=1=15 

120=1=14 

VCC1062GDL-10 

187.036740 

9.806686 

799=1=10 

100=t 7 

VCC1062GDL-11 

187.025115 

9.796717 

606=1=11 

138=1=10 

VCC1062GDL-12 

187.025085 

9.798357 

610=1=10 

116=1= 7 

VCC1062GDL-13 

187.000560 

9.796535 

386=1=10 

120=1= 9 

VCC1062GDL-14 

187.005885 

9.798541 

385=1=11 

101=t 8 

VCC1062GDL-15 

187.004535 

9.801900 

350=1=11 

llli 6 

VCC1062GDL-16 

187.001925 

9.804249 

346=1=10 

104=t 5 

VCC1062GDL-17 

187.003020 

9.802365 

333=1= 9 

102=1= 6 

VCC1062GDL-18 

186.999240 

9.805346 

356±12 


VCC1062GDL-19 

186.995685 

9.807220 

352=1=10 


VCC1062GDL-20 

186.994635 

9.807293 

362=1=10 


VCC1062GDL-21 

187.037685 

9.798342 

665±12 


VCC1062GDL-22 

186.992850 

9.808851 

382d=ll 


VCC1062GDL-23 

186.991335 

9.802511 

340±10 


VCC1062GDL-24 

186.986835 

9.806428 

381±12 


VCC1062GDL-25 

186.985650 

9.798775 

376=1=19 


VCC1062GDL-26 

187.012890 

9.792312 

501±19 


VCC1062GDL-27 

187.028310 

9.813976 

608±17 


VCC1062GDL-28 

187.018665 

9.791675 

527±25 


VCC1062GDL-29 

187.012260 

9.817571 

572±25 


VCC 685 

VCC685GDL-01 

185.985465 

16.692770 

1224±34 

129=1=27 

VCC685GDL-02 

185.988135 

16.683411 

1042=1=11 

99=1=12 

VCC685GDL-03 

185.997015 

16.701620 

1381d=12 

94=t 7 

VCC685GDL-04 

185.994150 

16.696199 

1372d=12 

103=1= 7 

VCC685GDL-05 

185.989155 

16.688009 

1034±12 

102=1= 9 

VCC685GDL-06 

185.991735 

16.694599 

1278±12 

136=1= 9 

VCC685GDL-07 

185.978040 

16.679291 

1060±90 


VCC685GDL-08 

185.989965 

16.701780 

1323=1=71 



Note. — Col.(l): Slit ID; Col.(2): Right ascension in decimal degrees (J2000); Col.(3): Declination in decimal degrees (J2000); Col.(4): 
Heliocentric radial velocity and uncertainty; Col.(5): Velocity dispersion and uncertainty. 













